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JUMPING AROUND THE CENTER OF K-SPACE 
IN ANNULAR SEQUENTIAL STEPS 

DESCailPTION 

The following relates to the magnetic resonance arts. It finds particular application 
in magnetic resonance imaging during a transient magnetic contrast agent bolus, and will 
be described with particular referrace thereto. However, it also finds application in 
magnetic resonance imagmg of transient events generally, such as in gated magnetic 
5 resonance imaging. 

In certain magnetic resonance imaging applications, k-space data is preferably 
acquired during a brief time interval. For example, in contrast-enhanced imaging the 
magnetic resonance imaging data are preferably acquired in the typically brief time interval 
during which the magnetic contrast agent bolus provides the most strongly enhanced 

10 magnetic resonance contrast in the tissue. In some contrast enhanced magnetic resonance 
angiography techniques, center of k-space data are collected while the contrast agent is m 
the arteries and before it arrives in the veins to provide good visualization of arterial flow 
by suppressing the not yet enhanced veneous signal. As another example, in cardiac gated 
imaging the magnetic resonance imaging data are preferably acquired in a short time 

15 interval corresponding to a selected cardiac phase so as to reduce motion artifacts. 

It will be appreciated that in these and similar applications the time interval for 
acquiring data is brief, but generally not well-defined. A contrast agent bolus, for example, 
typically produces a ramp-iq) in magnetic resonance signal intensity as the contrast agent 
washes into the imaged region. The slope of the signal intensity rise begins to level off 

20 until eventually a peak magnetic resonance signal intensity is reached Beyond this 
intensity peak, the signal intensity begins to decrease as the contrast agent washes out of 
the imaged region. Hence, the brief time interval for optimal imaging is not well-defined, 
but it is usually advantageous to concentrate acquisition of data around the time of the 
maximum bohis signal intensity. For example, separation of arterial and veneous signal is 

25 achieved by acquiring the center of k-space data while the bolus induced signal maYirr^iim 
occiurs in the vessel of interest. Similarly, in gated cardiac imaging, the objective is to 
image a single cardiac phase; as the imaging tune increases, the cardiac muscle 
increasmgly deviates from the target cardiac phase. 



i 
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To improve imaging of a transient event, it is known to order the k-space sampling 
so that k-space samples around the center of k-space are acquired closest to the optimal 
time, with k-space samples away from k-space center being acquured at less optimal times. 
For example, Riederer et al., U.S. patent no. 5,122,747 discloses a type of centric sampling 
in which the sampling trajectory spirals outward from k-space center, so that the initial 
k-space samples acquhed at the optimum time are near k-space centra-. The k-space 
samples close to the center of k-space are low frequency samples that contribute 
substantially to the overaU image. By distmction, k-space samples at larger distances from 
the center of k-space are high frequency samples that contribute to fino- details of the 
image. Hence, by acquiring san^les near the center of k-space first, the samples that most 
contribute to the overall image are acquired closer to the botes peak, target arterial or 
veneous phase, target cardiac phase, or other optimal time iniervaL 

Ahhough data acquisition techniques exist which advantageously concentrate data 
acquisition toward k-space center, these existing techniques have certain disadvantages. 
The data are not collected using a geometry that conqwrts with typical Fourier 
transform-based reconstruction algorithms. Hence, a substantial amount of sorting of 
k-space samples is usually required to prepare the collected imaging data for image 
reconstruction processmg. Such sorting can be computationally intensive, slowing down 
the data processing and delaying generation of useful images. 

Moreover, by collecting data spiraling out from the center, substantially tiie entire 
data set must be collected before tiiere is a completely acquired plane of k-space on which 
to begin Fourier-based image reconstruction. This prevents concurrent data acquisition and 
image reconstmction processing, which ftirther delays the generation of usefiil images. 

Yet another problem with existing techniques is that they typically involve 
applying a substantial number of large magnetic field gradient steps in order to follow the 
complex centric or other selected k-space acquisition trajectory. These large magnetic field 
gradients can generate eddy currents that contribute to image artifects. 

The present mvaition contemplates improved apparatuses and methods that 
overcome the aforementioned limitations and others. 
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According to one aspect, a magnetic resonance imaging method is provided. A 
k-space is divided into a central region disposed at k-space center and one or more annular 
surroimding regions having increasing distances from k-space center. The one or more 
5 annular surrounding regions including an outermost surrounding region having a largest 
distance from k-space center. The k-space saixq>les in the central region are acquired. 
Subsequent to the acquiring of k-space samples in the central region, k-space san^)les in 
the one or more annular surrounding regions (104, 106) are acquired. The k-space samples 
in the outermost surrounding region (106) are acquired last The k-space sanq>les are 

10 acquired in at least the outermost surrounding region using a row-by-row data acquisition 
ordering in which each row of k-space samples acquired in the outermost surrounding 
region, together with selected already-acquired k-space data from the regions oflier than the 
outermost surroimding region, forms a completed data set for reconstructing an image 
plane. Each completed data set is reconstructed into a reconstructed image plane without 

IS waiting for all k-space samples in .the outermost surrounding region to be acquired such 
that the reconstructing occurs at least partially concurrently with the acquiring. 

According to anoth^ aspect, a magnetic resonance imaging apparatus is disclosed. 
A magnetic resonance imaging scanner images an associated imaging subject A magnetic 
resonance imaging controller performing a method including: (i) dividing k-space into a 

20 central region disposed at k-space center and one or more annular surrounding regions 
having increasing distances from k-space center and including an outermost surrounding 
region of largest distance from k-space center; (ii) determining an optimum time for 
imaging a magnetic contrast agent bolus; (iii) acquiring k-space samples in tiie central 
region at about the optimum time, and (iv) after acquiring the k-space samples in the 

25 central region, acquiring k-space samples in the one or more annular surrounding regions. 
The acquiring in at least the outermost surroundmg region uses a plane-by-plane data 
acquisition ordering in which all k-space samples in the outermost surrounding region 
belonging to a current k-space plane are acquired to complete the current k-space plane 
before samples in the outermost surrounding region belonging to other k-space planes are 

30 acquired. A reconstruction processor reconstructs the completed cmrent k-space plane into 
a reconstructed plane image without waiting for other k-space planes to be completed. 
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According to yet another aspect, a magnetic resonance imaging apparatus is 
disclosed A means is provided for dividing k-space into a central region disposed at 
k-space center and one or more annular surrounding regions having increasing distances 
from k-q)ace center. The one or more annular surrounding regions include an outemiost 
surrounding region having a largest distance from k-space center. A means is provided for 
acquiring k-space samples in the k-space. The k-space samples in fhe central region are 
acquired first The k-space samples in the outermost surrounding region are acquired last. 
The k-space samples in at least the outermost surrounding region are acquired using a 
row-by-row data acquisition ordering in whidi each row of k-space samples acquired in 
the outennost surrounding region conqiletes a k-space plane. A means is provided for 
reconstructing each conq)leted k-space plane into a reconstructed image plane without 
waiting for all k-space san[q>les in the outennost surrounding region to be acquired. 

One advantage resides in shortened imagiTig session times. 

Another advantage resides in reduced image artifacts. 

Another advantage resides in enabling concurrent imaging data acquisition, 
reconstruction, and display processing. 

Yet another advantage resides in reduced data sorting complexity. 

Still yet another advantage resides in reduced generation of eddy ciurents. 

Numerous additional advantages and benefits will become apparent to those of 
ordinary skill in the art upon reading the following detailed descrq)tion of the preferred 
embodiments. . 

The invention may take form in various components and arrangements of 
coroponents, and in various process operations and arrangements of process opemtions. 
The drawings are only for the purpose of illustrating preferred embodunents and are not to 
be construed as limiting the invention. 
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FIGURE 1 diagraiDinatically shows a magnetic resonance imaging system. 

FIGURE 2 shows the division of the k-space sampling into (i) a central region 
disposed at k-space center, and (ii) two annular smroimding regions. 

FIGURE 3 shows a suitable random or pseudorandom acquisition ordering for 
5 k-space samples lying in the central k-space region. 

FIGURE 4 shows a suitable row-by-iow acquisition ordering for k-space samples 
lying in the inner annular surrounding k-space region. 

FIGURE 5 shows a suitable row-by-row acquisition ordering for k-space samples 
lying in the outermost annular k-space region. 
10 FIGURE 6 diagrammatically shows relative timing of the data acquisition, the data 

reconstruction, and the data display aspects of the imaging. 

FIGURE? shows a suitable embodiment of the reconstruction processor of 
nGURE 1. 

FIGURE 8 shows a suitable random or pseudorandom acquisition ordering for 
IS k-space samples lying in the inner annular surrounding k-space region. • ' 

With reference to FIGURE :1, a magnetic resonance imaging scanner 10 includes a 
housing 12 defining a generally cylindrical scanner bore 14 inside of which an associated 
imaging subject 16 is disposed. Main magnetic field coils 20 are disposed inside the 

< 

housing 12, and produce a main Bo magnetic field directed generally parallel with a central 
20 axis 22 of the scanner bore 14. The main magnetic field coils 20 are typically 
superconducting coils disposed inside cryoshrouding 24, although resistive main magaets 
can also be used. TThe housing 12 also houses or supports magnetic field gradient coils 30 
for selectively producing magnetic field gradients in the bore 14. The housing 12 further 
houses or supports a radio firequency body coil 32 for selectively exciting and/or detecting 
25 magnetic resonances. The housing 12 typically includes a cosmetic inner liner 36 defining 
the scanner bore 14. For some imaging applications, one or more local radio frequency 
coils such as a surface coil 44 may be used for radio frequency excitation and/or detection. 

In contrast-enhanced magnetic resonance imaging, a bolus injector 46 is used to 
rapidly infuse a selected magnetic contrast agent, such as a gadolinium-based contrast 
30 agent, into the imaging subject 16. For a bohis injection into the bloodstream of a human 
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imaging subject, the contrast agent rapidly washes into and out of tissue of interest. The 
bolus transit time between wash-in and wash-out is typically on the order of a few seconds. 
In the iUustrated embodiment, the bolus injector 46 is coupled into the bloodstream of the 
imaging subject 16 by a catheter 48. 
5 A magnetic resonance imaging controller 50 operates magnetic jQeld gradient 

controllers 52 to selectively energize the magnetic field gradient coils 30, and operates a 
radio frequency transmitter 54 coupled to Ae radio frequency coil 32 as shown, or coig)led 
to a local coil such as the surface coil 44, to selectively inject radio frequency excitation 
pulses into the imaging subject 16. By selectively operating the magnetic field gradient 

10 coils 30 and the radio frequmcy coil 32, magnetic resonance is generated and spatially 
encoded in at least a portion of a region of interest of the imaging subject 16. During 
imaging data acquisition, the magnetic resonance imagmg controller 50 operates a radio 
frequency receiver 56 coupled to the radio frequency sur&ce coil 44 as shown, or coupled 
to the whole body coil 32, to acquire magnetic resonance san^)les that are stored in a 

15 magnetic resonance data buffer 60. 

The acquired imaging data are reconstructed by a reconstruction processor 62 into 
an image representation, typically using a Fourier transform-based reconstruction 
algorithm. The reconstructed image or ioaages generated by the reconstruction processor 62 • 
are stored in an images memory 64, and can be displayed on a display device 66 of a user 

20 interfrice 68, stored in non-volatile memory, transmitted over a local intranet or the 
Internet, viewed, stored, manipulated, or so forth. The user interface 68 also preferably 
enables a radiobgist, technician, or other operator of the magnetic resonance imaging 
scanner 10 to communicate with the magnetic resonance imaging controller 50 to select, 
modify, and execute magnetic resonance imaging sequences. 

25 The described magnetic resonance imaging system is an illustrative example only. 

The magnetic resonance data acquisition techniques described herein can be used with 
substantially any type of magnetic resonance imaging scanner, including but not limited to 
horizontal bore scaimers, vertical bore scanners, open magnet scanners, and so forth. 
Moreover, the magnetic resonance data acquisition techniques described herein can be 

30 enployed in conjunction with arrays of receive coils for parallel imaging, SENSE imaging, 
and the like. 
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In contrast enhanced imaging, the bolus injector 46 administers a magnetic contrast 
agent bolus to the imaging subject 16. When.administered into the bloodstream, the bolus 
washes into the imaging region of interest, and then washes out, typically on the order of a 
few seconds. Accordingly, the imaging data acquisition is preferably timed to acquire 
5 imaging data at the peak of the magnetic contrast provided by the bohis in the imaging 
region of interest As will be described, the imaging data is acquired in a manner in which 
k-space samples in a region near the center of k-space are acquired first, followed by 
samples in regions at larger distances firom k-space center. In one approach for timing the 
imaging, the magnetic resonance imaging controller 50 sends a bolus triggering signal 70 

10 to the bolus injector 46. The bolus injector 46 administers tiie bolus via the catheter 48 
responsive to receiving the bohis triggering signal 70. The magnetic resonance imaging 
controller 50 initiates fhe magnetic resonance imaging data acquisition at a pre-selected 
time after sending the bohis triggering signal 70, where the pre-selected time corresponds 
to the time it takes for the administered bolus to wash into the imaging region of interest 

15 In anoth^ approach for timing the imaging data acquisition, the scanner 10 excites 

and acquires a projection or other magnetic resonance signal firom the region of interest • 
The magnetic resonance unaging controller 50 monitors the intensity of this signal. The 
bolus injector 46 adndnisters the magnetic contrast agent bolus. As the administered bolus 
washes into the imaging region of interest, the intensity of the projection or other magnetic 

20 resonance signal changes accordingly. The magnetic resonance imaging controller 50 
detects this change and begins image data acquisition. In yet other embodiments, such as in 
cardiac gated imaging, the magnetic resonance imaging controller 50 uses a suitable 
cardiac gating signal, for instance extracted fipom an electrocardiograph signal, to trigger 
data acquisition. Similarly, substantially any other physiological event can be detected and 

25 used to trigger initiation of magnetic resonance imaging of a volume of interest affected by 
the physiological event 

With continumg reference to FIGURE 1 and with fiirther to FIGURES 2-5, the 
magnetic resonance imaging controller 50 divides k-space 100 into a centraFregion 102 
disposed at k-space cent^ and one or more annular surroimding regions having increasing 

30 distances firom k-space center, namely two annular surrounding regions 104, 106 in the 
example embodiment illustrated in FIGURE 2. The k-space representations of FIGURES 
2-5 show a view perpendicular to a primary data acquisition coordinate designated "Cp", 
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and perpendicular to an orthogonal secondary data acquisition coordinate designated "cs". 
An unlabeled readout coordinate is oriented transverse to illustrated the primary Cp and 
secondary Cs coordinates, and hence parallel with the direction of view. In some 
embodiments, the Cp coordinate is a slice-select coordinate while the Cs coordinate is a 
5 phase-encode coordinate; however, the primary and secondary coordinates can have 
substantially arbitrary spatial orientations. For example, the phase-encode direction could 
be designated the primary coordinate Cp while the slice-select direction could be designated 
the secondary coordinate Cs, or one or both of primary Cp and secondary Cs coordinates may 
be inclined to the conventional sagittal, coronal, and axial anatomical coordinates. 

10 Moreover, while in the illustrated embodiment the primary Cp, secondary Cs, and readout 
data acquisition coordinates are mutually orthogonal, in some contemplated mibodiments 
a non-orthogonal data acquisition coordinates Systran may be miployed. 

In FIGURES 2-5, a rectangular grid of k-space samples 110 is illustrated In these 
FIGURES, each k-space sample or readout line is depicted as a dot The illustrated grid is 

15 10x10 for simplicity of illustration. In practice, grids of 128x128 to 1024x1024 are 
routine. In embodiments in which a three-dimensional volume is imaged, these k-space 
samples 110 each correspond to a readout line of k-space data that extends along the 
readout coordinate (that is, extends transverse to the primary Cp and secondary Cs 
coordinates, and thus appears as a dot from the perspective of FIGURES 2-5). Although a 

20 square 10x10 grid of samples is illustrated, the numb^ of samples in the Cp and Cs 
directions can be different. In such three-dimensional volume imaging, the central region 
102 will extend into the readout direction to define a generally cylindrical central region, 
while the annular regions 104, 106 will similarly extend into the readout direction to define 
three-dimensional annular volumes surrounding the cylindrical central region volume 102. 

25 The magnetic resonance imaging controller 50 acquires the central region 102 first 

using an acquisition ordering shown in FIGURE 3, and then acquires the inner annular 
surrounding region 104 using an acquisition ordering shown in FIGURE 4, and finally 
acquires the outermost annular surrounding region 106 using an acquisition ordering 
shown in FIGURE 5. Some preferred sampling orderings for use in acquiring data from 

30 each of these regions 102, 104, 106 is next described 

With particular reference to FIGURE 3, eadi k-space sample 110 in the central 
region 102 is labeled with its (n:dinal number in the acquisition sequence. The example 
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central region 102 encompasses twelve k-space samples labeled with ordinal numbers 1, 2, 
3, ... 12. In the illustrated embodiment, a substantially random or pseudorandom ordering 
of the k-space samples acquisition is used for acquiring the k-space samples in the central 
region 102. A random or pseudorandom ordering for acquiring samples of the central 
5 region 102 advantageously provides robustness with respect to mors in triggering the data 
acquisition. If the magnetic resonance data acquisition is triggered slightly early or slightly 
late such that the peak intensity contrast provided by the magiietic contrast agent bolus 
does not precisely coincide with acquisition of the central region 102, then this randcm or 
pseudorandom ordering reduces the sensitivity of the image quality to this triggering eoor. 

10 In this context, "random or pseudorandom'* is intended to encompass^ in addition to the 
ordinary meaning of those tarns, any pre-determined ordering selected to have statistkal 
characteristics similar to those of a random or pseudorandom ordering. 

A random or pseudorandom ordering typically has some cost in terms of iuCTeased 
magnetic field gradients employed to follow large jumps between the random sampling 

15 positions. However, because the central region 102 is generally relatively small and close 
to the center of k-space, the magnetic field gradients required to make the random Jumps 
within the central region 102 are typically acceptable. Alternatively, a centric or other 
ordering that collects the k-space sanq>les closest to the center of k-space can be used in 
the central region 102 instead of the illustrated random or pseudoiandom ordering. 

20 With particular reference to FIGURE 4, once all the k-space samples in the central 

region 102 are acquired, the magnetic resonance imaging controller 50 acquires k-space . 
samples in the inner surroimding aimular region 104 using a serpentine row-by-row 
acquisition ordering that is designed to acquire a row of readout samples having a givra Cp 
coordinate value as a group (that is, designed to acquire a row of readout samples along the 

25 Cs coordinate direction) while limiting the magnitudes of the applied magnetic field 
gradients. The acquisition order of the k-space samples is represented in FIGURE 4 where 
each acquired row of samples is indicated by an arrow labeled fay its ordinal numb^ of 
plane acquisition. In FIGURE 4, there are six rows labeled 1, 2, 3, 4, 5, 6. The direction of 
the arrow indicates the direction the row of samples is traversed during data acquisition. 

30 The first acquired sample is at the bottom of FIGURE 4 around the center position of the 
primary Cp coordinate. Starting at this k-space sample, acquisition of a row of saiiq)les 
traverses a positive or upward direction in Cs coordinate, and does not include acquisitbn 
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of those k-space samples along the row which lie inside the central region 102 and which 
thus have aheady been required. It will be appreciated that traversing this row of k-space 
samples involves gradient steps in the secondary Cs coordinate direction, without 
application of any gradients in the primary Cp direction. 

Once the first row of k-space samples are acquired, a gradient pulse in the primary 
Cp direction is applied to move to another row, which also lies near the center position of 
the primary Cp coordinate. This second row of k-space sanaples is traversed in the negative 
or downward direction, that is, the second row is traversed in the opposite direction fiom 
the first row. The seccmd row is also traversed using only secondary Cg coordinate gradient 
pulses, without application of any gradients in the primary Cp direction. 

Once the second row k-space samples is acquired, a gradient pulse in the primaiy Cp 
directwn is applied to move to the tlurd row, which lies furflier out fiom the center position 
of the primary Cp coordinate. The third row of samples is traversed in the positive or 
upward directfon. Once the third row of k-space samples is acquired, a gradient pulse in 
the primaiy Cp direction is applied to move to the fourth row, which is traveled in the 
negative, or downward durection. Once the fourfli row of k-space samples is acquired, a 
gradient pulse in the primary Cp direction is applied to move to the fifth row, which lies at 
an edge of the surrounding annular region 104. The fifth row of k-space samples is 
traversed in the positive or upward direction. Once the fifth row of k-space samples is 
acquired, a gradient pulse in the primary Cp direction is applied to move to the sixth and 
final row of k-space samples intersecting the surrounding region 104, which lies at the 
opposite edge of the suirounding annular region 104. The sixth row of k-space samples is 
traversed in the negative or downward directioa 

The described row-by-row acquisition pattern m the inner surrounding annular 
region 104 has certain advantages. In three-dimensional Fourier-encoded imaging, each 
k-space sample 110 corresponds to a readout line of k-space m the readout direction. 
Hence, a row of k-space samples corresponds to a plane of k-space. This plane of k-space 
can be transformed into a reconstructed image plane using two-dimensional Fourier-based 
reconstruction processmg, thus elimmating data processing delays caused by complex 
sorting operations. Moreover, the row-by-row acquisition ordering limits the magnitudes of 
changes in amplitude of the magnetic field gradients used m traversing the trajectory, 
which reduces eddy currents that can lead to hnage arti&cts. 
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With particular reference to FIGURE 5, once all the k-space samples in the central 
region 102 and in the inner annular sunounding region 104 are acquired, the magnetic 
resonance imaging controller SO acquires k-space samples in the outermost surrounding 
annular region 106 using the same type of row-by-row acquisition ordering used in 
S acquiring samples in the inner region 104. The acquisition of each row of k-space samples 
is represented in FIGURE 5 by an arrow labeled by its ordinal number of acquisition. 
There are ten row of k-space samples labeled 1, 2, 3, 10. The direction of the arrow 
indicates the directbn the row of k-space samples is traversed during the data acquisition. 
In similar fashion to the acquisition ordering used in the inner surrounding region 104, the 

10 row of the outermost surrounding region 106 are acquired such that the first and second 
rows having primary Cp coordinate values closest to the center of the outermost region 106 
are acquired first, and rows with increasing primary coordinate value are subsequently 
acquired, finishing with acquisition of the ninth and tenth rows of k-space sanq)les located 
at the opposite edges of the outermost annular region 106, respectively. The direction of 

15 acquisition of sanq)les along a row in the secondary Cs direction reverses at llie beginning 
of each new row to limit applied magnetic jSeld gradient magnitudes. Thus, rows 1, 3, S, 7, 

9 are traversed in Ihe positive or upward direction, while rows 2, 4, 6, 8, 10 are traversed in 

» 

the negative or downward direction. When acquiring the k-space data of the outermost 
surrounding region 106,. k-space samples along the row which lie in the central region 102 
20 or in the inner annular surrounding region 104, and which have thus already been acquired, 
are not re-acquired during the row-by-row san^ling of the outermost surroxmding region 
106. 

With continuing reference to FIGURES 1-5 and with further reference to FIGURE 
6, the . row-by-row acquisition of the outermost annular surrounding region 106 

25 advantageously enables image reconstructiaB to initiate before all k-space data is acquired, 
and allows image reconstruction to be poformed substantially concurrently with data 
acquisition once a fiill row of k-space samples (corresponding to a plane of k-space) has 
been acquired. As diagrammatically shown in FIGURE 6, the k-space samples in the 
central region 102 are acquired during a time interval 120, followed by acquisition of the 

30 k-space san^les in the inner annular sunounding region 104 during a subsequent time 
interval 122. At the end of the second acquisition time interval 122, no plane of k-space 
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has been fiilly sampled, since the outennost annular region 106 includes at least some 
k-space samples from every plane of k-space in the imaging volume. 

The row-by*row sampling of the outermost annular region 106 is performed during 
a time interval 124 following the time interval 122. Once the first row of k-space samples 
is acquired, all the data for a fiill plane of k-space are available, and reconstruction of that 
plane can be initiated by the reconstruction processor 62. 

By way of example, when the primary Cp direction is the sUce-select direction, the 
secondary Cs direction is the phase-encode direction, and the imlabeled direction 
perpendicular to and Os is the readout direction, acquisition of a row of k-space samples 
: corresponds to acquiring all phase-encode steps for a selected slice. As soon as all of the 
Jlc-^ace samples along the first row of the out^most region 106 are acquired by traversing 
the phase-encode steps along the row for diat selected slice (jimiping over those 
phase-encode steps aheady acquired during acquisition of regions 102, 104), the central 
slice can be reoonstructed. When all of the data along row 2 is collected, the adjoining slice 
can be reconstnicted, and so forth. Moreover, the same principle applies for other 
two-dimensional planes through three-dimensional k-space. 

Hence, a reconstruction time interval 130 initiates subsequent to completion of 
sampling of the first row of k-space samples in the outennost aimular region 106. The 
reconstruction may extend beyond the time interval 124 during which the samples in the 
. outermost annular region 106 are acquired (as indicated by the curved right-side "bieak" of 
the reconstruction time interval 130 in FIGURE 6); however, substantial data processing 
speedup is obtained through the partial overlap of the acquisition time interval 124 and the 
reconstruction time inteml 130. 

With brief retuniing reference to FIGURE 1, in one suitable embodiment, the 
magnetic resonance imaging controller 50 sends a reconstruction trigger signal 72 to the 
reconstruction processor 62 each time a row of k-space samples has been con[q>letely 
acquired. Receipt of the reconstruction trigger signal 72 causes the reconstruction 
processor 62 to begin reconstruction of the corresponding image plane. In another 
approach, the reconstruction processor 62 monitors the sampling data buffer 60 to 
determine when acquisition of a row of k-space samples has been completed. 

With returning refermce to FIGURE 6, the data processing does not involve a 
substantial amount of data sorting. The surrounding regions 104, 106 are acquired on a 
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row-by-row basis that con5)orts with two-dimensional Fourier reconstruction, so that 
substantial sorting is not required for those data. The data from the central region 102 is 
randomized or pseudorandomced and it must be sorted; however, this sorting is performed 
contemporaneous with the acquisition of the data, with each sample placed into a position 
5 in the data buffer corresponding to its k-space coordinates as the next sample is being 
acquired. 

Still further, since the reconstruction processor 62 reconstructs each plane as its 
k-space data acquisition is oon^leted during the acquisition interval 124, the image of that 
reconstructed plane can be displayed by the display device 66 in a displaying time int^al 
10 136 that at least partially ovodaps the acquisition time interval 124 and the reconstiuction 
time interval 130. Thus, the image begins to be displayed even while data acquisition is 
ongoing. 

With reference to FIGURE 7, a suitable embodiment of die reconstruction 
processor 62 is described. In some embodiments, the sorting of the data is determined prior 
IS to the acquisition and the expected acquisition order is provided to the reconstniction 

« 

processor 62 in an array of k-space points determined by the magnetic resonance imaging 
controller 50. Hence, no sorting' is performed during reconstruction beyond the initial 
placement of each data saiiq>le into a location in the data buffer corresponding to its 
position in k-space by means of a data sorter 150. With reference to FIGURE 6, this 

20 sorting by the reconstruction processor 62 during the time interval 132 is performed 
concurrently with the acquiring and reconstructing of the data. 

In three-dimensional imaging, each k-space sample 110 corresponds to a readout 
line in the readout direction which is orthogonal to the primary Cp and secondaiy Cs 
directions. Accordingly, a one-dimensional Foxirier transform 152 is applied to each 

25 k-space sample to recover the spatial content in the readout direction. 

Within the reconstruction processor 62, a slice checker 160 determines when a 
conc^lete row of k-space samples, corresponding to a complete data set for reconstracting 
an image plane, has been acquired. Alternately, a reconstruction trigger signal 72 (see 
FIGURE 1) could indicate when an image plane data set has been acquired. The convicted 

30 data set, which has akeady been sorted by the data sorter 150 and Fourier transformed in 
the readout direction by the Fourier transform processor 152, is processed by a second 
one-dimensional Fourier transform processor 162 operating in the secondary Cs coordinate 
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direction to recover the spatial content along the secondary Cs coordinate direction. The 
output of the second one-dimensional Fourier transform processor 162 is a reconstructed 
image of the just-conqjleted primary Cp plane. As acquisition of each row of k-space 
samples corresponding to an image plane is completed during the time interval 124, that 
S just-conQ)leted image plane data set is similarly reconstructed 

The dividmg up and sampling of k-space illustrated in FIGURES 2-5 is only an 
example. In some contemplated embodiments, the inner surrounding region 104 is omitted, 
Le., the data are divided between a randomly acquired central region and a peripheral 
region that is acquired row-by-row. In other conten[q)lated embodiments, the number of 

10 surrounding annular regions is greater than two. As shown in FIGURES 2 and S, the 
^ example outermost annular region 106 has a circular or oval inner boundary and a square 
or rectangular oxster boundary (in the Cp-Cs plane). Thus, it will be appreciated that the 
annular surrounding regions can have various cross-sectional shapes. 

The advantageous at least partial concurrency of data acquisition and reconstruction • 

15 . (and optionaUy also concurrmt displaying) is obtained oiily during acqdsition of k-q)a 
samples in the* outermost region 106, that is, only during the time interval 124. 
Accordingly, there are advantages to making the outeraiost region large so that it contains 
a large fraction of the k-space saniples compared with inner regions 102, 104, thus 
enabling a large overlap of data acquisition and image reconstruction. 

20 On the other hand, the row-by-row acquisition empbyed in the outermost region 

106 is not strongly biased toward acquiring data near the center of k-space, and so making 
. the outermost region 106 too large can increase temporal artifacts. To reduce temporal 
artifacts, the central region 102 should be made larger so that central k-space samples that 
contribute strongly to the overall image are acquired first 

25 As the central region 102 is enlarged, however, the sorting time prior to data 

acquisition is increased, causing an increase in acquisition planning time. Moreover, if the 
central region 102 is too large it will contain k-space samples relatively distant from 
k-space center. Including one or more iimer surroimding annular regions such as the inner 
surrounding annular region 104, sampled on a row-by-row basis as shown in FIGURE 4, 

30 provides a balance by biasing the early portion of the acquisition toward k-space center 
while providing a row-by-row acquisition ordering that is conducive to rapid data 
processing. 
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In some embodiments, the size of the central region 102 is determined as follows. 
The magnetic resonance imaging parameters for the acquiring of k-space samples is 
determined These parameters include a data acquisition rate. The time interval 120 for 
acquiring samples in the central region 102 is also selected, for example based on an 
5 expected duration of the magnetic contrast agent bolus in the tissue of interest Based on 
the data acquisition rate and the selected time interval 120, the number of k-space samples 
of the central region 102 is determined, and a suitable central region 102 corresponding to 
this number of samples is selected 

With reference to FIGURE 8, in some embodiments it is contemplated to use an 

10 ordering other than serpratine row-by-row ordering for the surrounding annular regions 
other than flie outermost surrounding aimular region. FIGURE 8, for example, shows an 
example of a random or pseudorandom sampling ordering for the iimer surrounding 
annular region 104. Substituting the random or pseudorandom sampling of FIGURE 8 for 
the row-by-row ordering of FIGURE 4 can be advantageous for certain imaging tasks. For 

IS example, this approach provides additional robustness against triggering time euors by 
extending the range of the random or pseudorandom sampling further out away fix)m 
k-space center, while the separate random or pseudorai^dom sampling of each of the central 
region 102 and the inner surrounding region 104 ensures that the random or pseudorandom 
sampling covers all the k-space samples within the boundary of the central region 102 

20 before sampling outside of central region 102. However, it will be recognized that the 
random or pseudorandom sampling of the inner surrounding region 104 increases the 
required data sorting time. 

Although not illustrated, in some embodiments it is contemplated to employ a 
row-by-row acquisition ordering for the central region 102 and for the one or more 

25 surrounding annular regions 104, 106. Row-by-row ordering does not by itself strongly 
bias the early sampling toward k-space center; however, by using a relatively small central 
region 102 and one or more relatively small inner surrounding annular regions 104 biasing 
of the early sampling toward k-space center can be achieved. Moreover, since in these 
embodiments all regions including the central region 102 are sampled on a row-by-row 

30 basis, sorting delays in the subsequent data processing are substantially reduced. 

The invention has been described with reference to the preferred embodiments. 
Obviously, modifications and alterations will occur to others upon reading and 
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understanding the preceding detailed description. It is intended that the invention be 
construed as including all such modifications and alterations insofar as they come within 
the scope of the appended claims or the equivalents thereof. 



